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Abstract. Recent experimental studies of the antiproton-deuteron system at low energies have shown that 
the imaginary part of the antiproton-deuteron scattering length is smaller than the antiproton-proton one. 
Two- and three-body systems with strong annihilation are investigated and a mechanism explaining this 
unexpected relation between the imaginary parts of the scattering lengths is proposed. 



PACS. 25.43.+t Antiproton-induced reactions 

1 Introduction 

Recent measurements of antiproton-deuterium and anti- 
proton-helium annihilation cross-section at low antiproton 
momentum (p lab down to 35 MeV/c) Q| and of the shift 
and the width of the antiproton-deuteron atomic ground 
' state have provided the first information on interaction 
, of slow antiprotons with light nuclei. Both experiments 
gave unexpected but coherent results. The width of the IS 
state of the pD atom appears to be approximately equal 
to that for protonium. In the region of low antiproton 
momentum, the annihilation cross-sections in the pp, pD 
and p 4 He systems are approximately equal. 

The phenomenological analysis of these data || allows 
the imaginary parts of the S-wave scattering lengths for 
, the three nuclei to be extracted. The values was found to 
decrease with atomic number 

Im a sc (pp) = -[0.69 ± O.Ol(stat) ± 0.03(sys)] fm, 
Im a sc (pD) = -[0.62 ± 0.02(stat) ± 0.05(sys)] fm, 
Im a sc (p 4 He) = -[0.36 ± 0.03(stat)±°;^(sys)] fm. 

This result is in direct contradiction with the naive geo- 
metrical picture of annihilation which would suggest that 
the imaginary part of the scattering length (or correspond- 
ing width and cross-section) increases with the size of the 
nucleus. 

Only a few calculations of the atomic pD system ||, 
|5|,[| have been carried out prior to the measurements. In 
the first one Q , the authors used an approximate iterative 
method to obtain the optical pD potential starting from 
elementary NN and NN interactions. This method was ap- 
plied previously to atomic data for heavier nuclei (jjg. In 
the second calculation ||, the pD system interaction was 
studied by solution of the three body equations. In both 
approaches, the elementary NN and NN interactions was 
chosen in separable form (this limits the conclusions due 



to particular properties of the separable interaction). In 
the two articles, the same numerical result was obtained: 
the width of the IS pD atomic level appeared to be smaller 
than that for the pp system and the necessity to push the 
calculations beyond the simple impulse approximation was 
emphasized. The same conclusion was obtained in Q . In a 
third the article fl, the authors have investigated mostly 
different approximate solutions of the three-body prob- 
lem. Due to the strong annihilation the results was found 
to be not very sensitive to the three-body dynamics. How- 
ever, the shifts and the widths of the IS pD atomic level 
obtained in this article are quite different from the ones 
observed experimentally 

The aim of the present work is to understand why the 
simple geometrical picture does not work and to propose 
a mechanism, at least for pD system, to explain the un- 
expected behavior of the scattering length. 

Our explanation is based on two points. The first one 
is the fact that the imaginary part of the scattering length 
for the NN optical potentials with strong annihilation is 
mainly determined by the diffuscncss of the exponential 
tail of the potential, but not by the details of inner part of 
interaction. Therefore one can expect that Im a sc (pp) and 
Im a sc (pD) can be, at least, close to each other, in spite 
of the different structure of interactions in pp and in pD 
systems. This expectation is confirmed by our three-body 
calculation of Im a sc (pD). The second one is the increase 
of Im a sc (pp) relative to Im a sc (pD) due to a resonance 
in NN system. 

We shall be concerned with general features of anni- 
hilation in the pp and pD systems which are essential for 
explaining the observed phenomenon. There are two pos- 
sible methods for taking the annihilation into account: 
either to introduce the imaginary part to the potential 
(optical model) or to introduce explicitly effective annihi- 
lation channels (coupled channel model). The first possi- 
bility is more simple but leads to some problems related 
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to the non unitarity of this approach (see JTcj] and refer- 
ences therein). The second takes the unitarity into account 
explicitly but the numerical treatment becomes quite in- 
volved. Thus, to simplify the calculations only the optical 
model approach is considered. Furthermore, a very simple 
antinucleon-nucleon interaction is used. 

To simplify the numerical treatment we will calculate 
the scattering length a s corresponding to a pure strong 
interaction. The scattering lengths extracted from the ex- 
perimental data 

G>sc £1X6 modified by the Coulomb interac- 
tion. For the pp system, there is a simple phenomenolog- 
ical relation between a s , a sc , and the Bohr radius B: 

A = £ +c 

with complex constant C = 7.17 + zl.46 which is valid 
within 1% accuracy as was demonstrated in p~Tfl ■ If one 
supposes that this relation is the same for the pD system 
one obtains the following imaginary parts of the scattering 
lengths: 

Im a s (pp) = -(0.86 ± 0.04) fm, 
Im a s (pD) = -(0.80 ±0.06) fm. 

The real parts of the scattering lengths are taken from 
H , all errors are added quadratically. One can see that in 
absence of Coulomb forces the pD scattering length is still 
smaller than that for the pp system. 

The article is organized as follows. Section || presents 
the properties of simple two-body interactions in the case 
of strong annihilation. In Section |l the three-body pD 
problem is solved numerically by means of Faddeev equa- 
tions and a mechanism explaining the observed effect is 
found. Finally, we provide a brief summary of the results. 



2 Two-body interactions 

Before calculating the pD scattering length obtained in the 
next section by the numerical solution of the three-body 
problem it is instructive to present some analytical results 
on the scattering length in a few simple two-body optical 
potentials. These examples emphasize that the diffuseness 
of the interaction determines the imaginary part of the 
scattering length. This remains valid for pD interaction. 
For the square well potential 



U(r) 



-U , r < R 
0, r>R 



the scattering length is equal to 

tan pi? 



R 1 



pR 



(1) 



Here p 2 = 2mUo with Uq = We t(p . In the case of pure 
imaginary optical potential, ip = ir/2. Hereafter, we sup- 
pose that ip > 0. Potentials with an attractive real part 
have < tp < tt/2, and repulsive potentials have 7r/2 < 



ip < tt. For a pure imaginary potential, the imaginary 
part of the scattering length as a function of the imagi- 
nary depth starts from zero, reaches a maximum and, in 
the limit of very strong annihilation (W — > oo), tends to 
zero as \ j\fW . The scattering length becomes real. This 
result is well-known in physics of hadronic atoms |fl2f . The 
scattering on this potential is equivalent to scattering on 
hard sphere. In this limit, the wave function tends to zero 
within the region of annihilation. This is the so-called "S- 
wave suppression" which appears in the optical potential 
with strong annihilation. 

Even in the limit when the range of interaction be- 
comes very large (R — > oo), the imaginary part of the 
scattering length does not go to infinity: 



lim Im a s 

R^oo 



1 V 
; COS — . 

V2mW 2 



The parameter which defines the strength of annihilation 
is a product \pR\ (see ([!])). So hereafter, the strong anni- 
hilation limit should be taken as implying \pR\ S> 1. 

The fact that the imaginary part of the scattering 
length becomes real when the annihilation is strong is re- 
lated to the sharp edge of the square well potential. For 
another analytically solvable potential 

U(r) = -U a c- r/r ° 

with Uq = We i>p , the scattering length is equal to 



2r 



7 + In z 



ttN (2z ) 



2 J (2z ) 



Here 7 = 0.577... is Euler constant, z 2 = 2m[/o?~Q, Jo(z) 
and Nq(z) are Bessel and Neumann functions of index 0. 

In the limit of strong annihilation (|zo| 3> 1), Im cl s is 
defined by the diffuseness of the interaction vq: 

lim Im a s = — ro(7r — ip). 
\zo\—>oo 

The same asymptotic result for Im a s is valid for any 
complex potential which has an "exponential tail" , for in- 
stance, a Woods-Saxon potential which can be considered 
as a first approximation to describe the antiproton-nucleus 
atomic data [|l3| . More sophisticated potentials containing 
more complicated dependencies are also discussed in the 
literature (see JT^J and references therein). 

In the limit of strong annihilation, Im a s depends nei- 
ther on the range of interaction R nor on its particular 
form within this region, but it is determined only by the 
diffuseness of the exponential tail r$. The wave- function 
if - exponentially decreases when r crosses the region of 
diffuseness of the potential: 



y/2m\U(p)\dp 



~ exp 



From the expression for the scattering amplitude: 



/(*) 



2m 



. (kr) U(r)<F(r)dr 
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(k being the center of mass momentum) one can see that 
the contribution of the inner part of the interaction is ex- 
ponentially small due to the behavior of the wave- function. 
Qualitatively this means that the particle in a strong imag- 
inary potential "annihilates" at a large distance and that 
the probability to penetrate to the origin is exponentially 
small. The scattering takes place on the exponential tail 
and is thus determined by rQ. 

As was shown in jl5[] , this phenomenon leads to the 
effect of saturation of the antiproton-nucleus and kaon- 
nucleus atomic widths. 

Another prediction can be made for the value of the 
antinucleon- nucleus annihilation cross-sections. At very 
low energies, where the annihilation cross-section a a is 
determined by Im a s : ka a « 4ir |Im a s \, its value does not 
depend on the atomic number (for nuclei with the same 
diffuseness) . Strictly speaking, this simple conclusion is 
true for n only (for p, the low energy expression for the 
annihilation cross-section depends on the real part of the 
scattering length and contains higher partial wave param- 
eters [HQ, so that this case needs more careful consider- 
ation) . 

To complete this analysis let us consider an elementary 
interaction in a separable form V = g\a)(a\ with complex 
coupling constant g (as chosen in The imaginary 

part of the amplitude as a function of g has the same 
qualitative one-bump behavior as the imaginary square 
well potential: 

5 |a(k = 0)| 2 



The scattering length for this potential is given by 



l + 2m.g/d 3 k|a(k)| 2 /fc 2 ' 

where a(k) = (k|a). 

For finite complex g, a s is complex also. However, for 
a very strong interaction \g\ — > oo, the scattering length 
becomes real for any form- factor (k|a). The reason for 
such behavior is the "poor" spectral properties of sepa- 
rable interactions: a rank-one separable potential is able 
to produce only one bound state and all dynamics is gov- 
erned by the position of the corresponding poles of the 
S-matrix. 

Till now, only potentials which have the same range of 
the real and of the imaginary part were considered. How- 
ever, the optical potentials describing elementary anti- 
proton-proton interaction are quite different. For instance, 
in the case of the popular Kohno-Weise potential jl?], 
|H| which reproduces rather well all data on low energy 
nuclcon- antinucleon interaction, the range of the real part 
of the interaction is larger than that of the imaginary part 
(see Figure 2 in Q). For these potentials, the existence 
of an external real part of the potential can significantly 
modify the situation. 

As an example, let's consider the simplest potential 
composed of two square wells where the internal potential 
is purely imaginary and external one real: 

-iW, r < i?i 
U(r) = { -V, Rx < r < R 2 . 
0, r > R 2 



1 

P 



tan 



p(R 2 — Ri) — arctan 



p tan kRi 



Here k 2 = 2imW and p 2 — 2mV. 

If the external part of the potential is strong enough 
to create nearthreshold bound or resonant states, one ob- 
tains: 



-Im a s 



1 



1 



^/2mW cos 2 p(R 2 -Ri) 



These states manifest themselves as an additional enhance- 
ment factor l/cos 2 p(i?2 — Ri)- 

In the case of an exact resonance (the bound state has 
an energy exactly equal to zero p(R 2 — Ri) = (2n + 1)^) 
instead of a decreasing function seen previously we obtain 
an increasing one when W — > oo 



-Im a. 



\fmW 



The presence of near threshold poles (bound states or res- 
onances) can significantly change the behaviour of Im a s . 



3 pD system 

In the case of pD scattering, it is possible to obtain numeri- 
cal solution of the Faddeev equations. In this approach, all 
three body effects are taken into account. We performed 
such a calculation with a simple model of the NN inter- 
action to study the effect of strong annihilation and find 
out whether it is possible to reproduce the experimentally 
observed tendency. All interactions in this three body sys- 
tem were chosen in local form to avoid problems related 
to particular properties of separable interaction. 

We solve the Faddeev equations in the coordinate space. 
The method can be found in pSi^ . Here we mention only 
main points. 



3.1 Description of the model 

The three body wave function in the Faddeev approach is 
represented as the sum of three components 

f = <fi(xi,yi) +<f r 2 (x 2 ,y 2 ) +#3(x 3 ,y 3 ). 

where x^, y^ {i = 1,2,3) denote three sets of Jacobi coor- 
dinates defined according to 

X ? ; Vj Vfc , 



2 {rj+r k 

Here is nucleon (or antinucleon) position vector. 
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The three-body Schrodinger equation is then rewritten 
as the system of equations for the Faddeev components: 

[k 2 + A X1 + A y 1 - mVi (xi )] W i (xi , y r ) = 
mVi (xi ) [V 2 (x 2 , y 2 ) + #3 (x 3 , y 3 )] , 

[k 2 + A X2 + A y2 - to^ 2 (x 2 )] ^ 2 (x 2 ,y 2 ) = 
mV 2 (x 2 )[^ 3 (x 3 ,y 3 ) +^i(xi,yi)], 



[k 2 + A X3 + A y3 - ™F 3 (x 3 )] ^3(x 3 ,y 3 ) - 

mV 3 (x 3 ) ( Xl , yi ) + V 2 (x 2 , y 2 )] . 



where E is the 3-body energy (above the mass of the 
Nnp system) and k 2 = raE. To simplify the treatment, 
we choose the NN interaction in the isospin independent 
form. Therefore we put V\(x) — V 2 (a;) = U(x) where U 
is the complex antinucleon-nucleon potential. In this case 
the components <?i and $2 have the same functional de- 
pendence on their arguments. We put also V^{x) = V d (x) 
where V d {x) is the np potential responsible for describing 
the deuteron. For V d (x) we use the Mafliet-Tjon poten- 
tial ||c]] . It corresponds to the S-wave np interaction only. 
Since we consider the low energy scattering, we keep also 
only the S-wave in the antinucleon-deuteron scattering. 
Hence, the Faddeev components depend on the moduli 
Xi,yi. Instead of &i(xi, yi) we introduce il>i(xi, yi): 



&i{xi,yi) 



tpi{xi,yi) 



•£%yi 



In the r.h. -sides of the equations (g), the moduli of the 
Jacobi coordinates x 2 , y 2 and £ 3 , y 3 expressed in terms of 
X! , y! , depend on the angle 8 between X! and yi . There- 
fore the partial wave decomposition results in the integrals 
over u = cos 9. In this way, the equations (Q) turn into the 
two component system of integro-diffcrcntial equations: 

dx 2 

— 2 - mV d {x) 



mU(x)~ 



8?_ 

dy 2 
1 



mU(x) 



dx 2 dy 2 



mV d (x) 



ipz{x,y) = 



du——ipi(x',y'), 
x y 



(3) 



with 



y 



2V3xyu + 3y 2 ) 1/2 , 
- 2VSxyu + y 2 ) 1/2 . 



We have taken into account that the dependence of the 
Faddeev components in the integrands on two sets of vari- 
ables (e.g., on X2, y 2 and x 3 , y 3 in the first equation in (||)) 
is reduced, because of integration over u, to the single set 
•<•'•. </'• 

Below the three-body threshold, in the asymptotic re- 
gion of x,y the component ip\{x,y) decreases exponen- 
tially. The component ^(x.y) describes the system asymp- 
totically (when the distance V3y/2 between antiproton 



and deuteron is large): 

i>l{x, y -> oo) = cp d (x) [sin(yq) + f(q) exp(iqy)} . (4) 

Here tp d (x) is the deuteron wave function, f(q) is the 
antiproton-deuteron scattering amplitude which is related 
to the corresponding phase shift 8 by f{q) = (exp(2iS) — 
l)/(2i). The momentum q is related to the center-of-mass 

pD energy as E%m — q 2 /m. With these definitions of /(g) 



(2) and q the scattering length is found to be: 



V3.. f(q) 
a = hm . 

2 g^o q 

The factor y/3/2 appears since it enters the relation be- 
tween the Jacobi coordinate y and the relative pD dis- 
tance. 

The NN potential was taken in the Woods-Saxon form 
to reproduce main features of Kohno-Wiese optical poten- 
tial 0,0: 



U(r) 



-V 
-iW 



1 + exp 
1 + exp 



r-R r 
r r 

r-R 

r,: 



with Rr = 1.2 fm, r r = 0.3 fm, R, = 0.55 fm, r t = 0.2 
fm (to reproduce the fact that the range of the real part 
is bigger than that of the imaginary part). V and W are 
considered as parameters. The imaginary part of the in- 
teraction is exactly the same as in the original potential 
0@- 



3.2 Pure annihilation and resonance limits 

The simplest and very instructive case of the NN interac- 
tion is the limit of pure annihilation V = 0. The imaginary 
part of the scattering length as a function of the imaginary 
depth W in two (NN) and three (ND) body problems are 
presented in Figure 

The results coincide with naive geometrical estima- 
tions: the imaginary part of the scattering length appears 
to be higher in the three body problem. For strong an- 
nihilation, the difference between the two values becomes 
smaller. This result is quite natural when we take into 
account the arguments given in the Section 2: for large 
values of W, Im a s does not depend on the size of nucleus 
R and is determined only by diffuseness. 

This limit shows that the observed relation between 
Im a s for pp and pD interactions cannot be explained 
within a model of pure annihilation. 

To make Im a s for the pp system bigger than for the 
pD system, it is necessary to create a near threshold res- 
onance or bound state (a pole of the S-matrix) in the two 
body system. Figure ^ shows an example when the real 
part of the two body interaction V is strong enough (80 
MeV) to be able to produce a S-wave pole which is very 
close to the threshold. Without annihilation (W = 0), one 
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0.5 - 



3 body problem 



2 body problem 



Fig. 1. The imaginary part of the scattering length as a func- 
tion of the imaginary part W of the Woods-Saxon potential 
when two body NN interaction is pure annihilation one (two 
body problem - dashed line; three body problem - solid line) 



has a loosely bound state with the S-matrix pole posi- 
tion in fc-plane of fcf, = i31.7 MeV/c. When annihilation is 
switched on this state acquires a width and is shifted away 
from the threshold: for W — 50 MeV, k b = -17.3 + i32.2 
MeV/c and for W = 500 MeV, k b = -130.9-M72.9 MeV/c. 

This state manifests itself by quite a large value of 
Im a s , as we discussed previously, in both two and three 
body problems. However, this resonance is less pronounced 
in the three body case and thus Im a s in this case is less 
sensitive to variation of the annihilation strength. 
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Fig. 2. The imaginary part of the scattering length as a func- 
tion of the imaginary part W of the Woods-Saxon potential 
when two body NN resonances are close to the threshold 



This example shows how one can change the relation 
between Im a s in two and three body problems. 

There are, at least, two ways to obtain a more realis- 
tic value of the imaginary part of the scattering length. 
Firstly, one can shift this state quite far from the thresh- 



old. Secondly, one can create not ground but an excited 
state which is sufficiently less sensitive to annihilation be- 
cause of the zero in the wave function. Figure |^ shows 
an example of such a situation (V = 500 MeV). With- 
out annihilation, the potential creates a loosely bound 
near threshold state (with k b = i36.5 MeV/c) which is, 
in fact, the first excited state (the ground state is very 
deep kb = i485 MeV/c). When the annihilation is switched 
on, the corresponding pole is shifted out of the threshold 
as previously but it is less sensitive to annihilation: for 
W = 50 MeV, k b = -11 + i37.8 MeV/c and for W = 500 
MeV, k b = -46 + i80.3 MeV/c. 
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Fig. 3. The imaginary part of the scattering length as a func- 
tion of the imaginary part W of the Woods-Saxon potential 
when two body NN resonances are far from the threshold 

It is interesting to note that in Kohno-Weise model [|l7| 
[Lsf such quite large resonances appear in all 4 (!) partial 
waves (with two different spins and isospins) pQ| . More- 
over, one of them is an excited state. 

As was found in |2^| these states appear also in a sep- 
arable potential which was used for the three body calcu- 
lations 



4 Conclusions 

Two- and three-body systems with very strong annihila- 
tion have been investigated using the optical model po- 
tential to understand the behaviour of the imaginary part 
of the scattering length as a function of the annihilation 
strength. 

In two-body systems, the imaginary part of the scat- 
tering length is sensitive mostly to the diffuseness of the 
potential. This circumstance gives a quite strong experi- 
mental prediction: the low energy antineutron-nucleus an- 
nihilation cross-section should be approximately the same 
for all nuclei. 

Due to Coulomb forces, the situation in antiproton- 
nucleus annihilation can be more complicated: higher par- 
tial waves and the real part of the scattering amplitude 



6 



V.A. Karmanov, K.V. Protasov, A.Yu. Voronin: Antiproton-deuteron annihilation at low energies 



could modify this simple picture. However, the experi- 
ments on the future AD (antiproton decelerator) facility 
would be of vital importance to understand the antinucleon- 
nucleus interaction. 

For three body systems, we propose a quite simple 
mechanism which is able to explain the observation that 
the imaginary part of the antiproton-deuteron scattering 
length is smaller than the antiproton-proton one. This 
mechanism is based on two points. Firstly, the imaginary 
part of the scattering length in both systems are deter- 
mined by the diffuseness of the annihilation potential. 
Thus these_ values are close to each other. Secondly, the 
two-body NN interaction produces S-wave poles and en- 
hances the values of the scattering length. This enhance- 
ment is more important in the pp system than in the 
pD one. These S-wave poles (resonances) are necessar- 
ily present in all models which describe the experimental 
data, however, their experimental observation is a delicate 
task because of their large widths. 
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